Universal phase boundary shifts for corner wetting and filhng 
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The phase boundaries for corner wetting (fiUing) in square 
and diagonal lattice Ising models are exactly determined and 
show a universal shift relative to wetting near the bulk criti- 
cality. More generally, scaling theory predicts that the filling 
phase boundary shift for wedges and cones is determined by a 
universal scaling function Rd{tp) depending only on the open- 
ing angle 2^p. Rd{tp) is determined exactly in d = 2 and 
approximately in higher dimensions using non-classical local 
functional and mean-field theory. Detailed numerical transfer 
matrix studies of the the magnetisation profile in finite-size 
Ising squares support the conjectured connection between fill- 
ing and the strong-fiuctuation regime of wetting. 



PACS numbers: 



;.45.Gd, 68.35.Rh, 68.45.- 



The wetting of structured and chemically heteroge- 
neous substrates has recently drawn a great deal of inter- 
est ||l|-|lO|] . This is motivated not only by the possible rele- 
vance to emerging technologies such as micro- fluidics 0, 
but also from a more fundamental statistical physics per- 
spective since surface structure may induce new types of 
interfacial phase transition. Examples of these are mor- 
phological phenomena on heterogeneous substrates [^|J^] 
and filling transitions for fluids adsorbed near corners, 
wedges and cones |]^-|l0||. Recent continuum effective 
interfacial Hamiltonian studies indicate that the condi- 
tions for observing continuous wedge and cone filling are 
much more relaxed than for continuous (critical) wet- 
ting. These studies also indicate that fluctuation effects 
at filling are characterised by a high degree of universality 
which persists even in the presence of long-ranged forces. 
For example, the filling of a three dimensional wedge 
is characterised by a universal roughness exponent ||] 
whilst for two dimensional corners (in both ordered and 
disordered systems iQ,^) and three dimensional cones pC| ] 
fluctuation effects mimic, quite precisely, behaviour pre- 
dicted to occur for the strong-fluctuation (SFL) regime 
of critical wetting ]T^ . 

Here we go beyond the effective Hamiltonian approach 
and investigate surface phase diagrams and fluctuation 
effects at filling using microscopic theory. We present ex- 
act results for two dimensional {d = 2) Ising models 
which improve upon previous low temperature, solid-on- 
solid (SOS) studies 0,11^ and show that, near the bulk 
critical point, the shift of the phase boundary (relative 
to wetting) is universal, determined only by the corner 
opening angle. The results of detailed numerical transfer- 
matrix studies, which test the proposed connection be- 



tween two dimensional filling and the SFL, are discussed 
at length and predictions for phase boundary shifts for 
wedges and cones in d — 3 are made. 

Consider a two dimensional corner with opening an- 
gle 2ip (with ip e (0,7r/2]) in contact with a bulk phase 
A which is at two-phase bulk AB coexistence. Thus 
ip = 7r/4 is a right angle corner whilst ip = 7r/2 is a 
flat, semi-infinite interface. Thermodynamic arguments 
|l6| indicate that the corner is completely filled by phase 
B when the contact angle 6 satisfies 



0{T) 



(1) 



implying that filling precedes wetting. The same con- 
dition applies for the filling of three-dimensional wedges 
and cones and is confirmed by effective Hamiltonian stud- 
ies 1^-^- Now recall that, for systems with sufficiently 
short-ranged forces, the scaling of the excess surface- free 
energy near the bulk critical temperature Tc [ p7[ implies 
that 9 = W{ahit-^^), where t = {T^ - T)/Tc is the re- 
duced temperature, hi is the surface field and Ai the 
surface gap exponent. Note that W{x) can be identi- 
fied as a universal scaling function provided we allow 
for a non-universal metric factor a. It follows that near 
the bulk critical point the filling phase boundary satisfies 



ah({T)t-^' 



W ^(f — "0) generalising the well known 



result for the wetting phase boundary (denoted hi{T)) 
corresponding to -tp = t:/2 ||T^. Consequently the phase 
boundary shift defined as 



T^Te hf{T) 



(2) 



should be universal depending only on the opening an- 
gle and dimensionality. Implicit in this definition is the 
assumption that both the wetting and filling transitions 
are continuous. This requirement is not restrictive in the 
two dimensional systems, which are the main concern of 
this paper, since both the filling and wetting transitions 
are continuous at all temperatures. In three dimensions 
we emphasise that both wedge and cone filling may be 
continuous even if the wetting transition is first-order. 
To continue, it follows that for fixed, small values of the 
surface field 



Tc - Tf 



l/Ai 



(3) 



which provides a useful means of estimating Tf provided 
Tu, is known. Moreover, whilst the character of the filling 



I 



transitions for the d = 3 wedge and cone are different, 
the focation of the phase boundary obeys the same ther- 
modynamic condition (|^) and the shift function R3{ip) is 
unique. 

To test the scahng theory we consider square and di- 
agonal lattice right-angle Ising corners which are natural 
extensions of Abraham's semi-infinite model |]l8| , |l9| (see 
Fig.l). We assume isotropic interactions in order to pre- 
serve reflection symmetry about the apex diagonal so 
that both boundaries behave as identical walls. The 
reduced interaction strength between nearest neighbour 
spins is K{= 1/T), except they are modified to 
hi/T) in the bottom row and far left column. The bound- 
ary spins are fixed to -1-1 whilst the bulk magnetisation 
(infinitely far from the boundaries) is supposed negative. 
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FIG. 1. Square (S) and diagonal (D) lattice isotropic Ising 
corners with weakened bonds. Typical low-temperature SOS 
interfacial configurations are shown. 

These boundary conditions induce a long domain wall 
which rounds the corner at some distance L (measured 
along the bottom row or left column) from the apex. On 
approaching the filling transition the mean value of L 
and hence the diagonal interfacial height z, measured in 
the [11] direction from the corner apex, diverge. Defin- 
ing the incremental free-energy Tp for each semi-infinite 
lattice as the difference of the wall-down spin and wall- 
up spin surface tensions recall that for square (S) and 
diagonal (D) Ising lattices Hpl 



cosh Tp = coshr 



1 



-F{T) 



(4) 



where = 2K + IntanhX and ^ 2 ln(sinh 2i4:) de- 
note the free interfacial tensions. The singular contribu- 
tions to the free energy are 



(e^(cosh2i^ - cosh2/\) - er*^ sinh2A') 
(cosh 2K - cosh 2k) sinh 2K 



(sinh^ 2K -1-2 sinh^ 2Kf 
sinh^ 2K 



(5) 



which vanish at wetting. Next note that for large L 
the excess free-energy of the corner must contain a term 
f{T)L arising from surface tension contributions with 



p 



T^/2 and /^(T) ^ tE' - r^. At filling 



f{Tf) = 0, leading to the exact conditions 

cosh 2k = cosh 2K - e"^^ sinh^ 2K 



sinh^2X^(^^^^i^^^^iKi^ 



-2K\ 



(6) 



(7) 



for the square and diagonal lattices respectively (see Fig. 
2). Both types of corner show the same universal phase 
boundary shift with 



i?2(7r/4) 



1 



V2+72 



(8) 



in agreement with scaling theory. The low temperature 
behaviour is also revealing. For the square lattice there 
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FIG. 2. Surface phase diagrams for the filling in the square 
and diagonal (shown on the left and right respectively) lat- 
tice Ising corners compared with critical wetting. The circles 
on the left hand figure show the numerically determined loci 
where the magnetisation profile shows finite-size scaling be- 
haviour characteristic of the SFL regime of wetting 

is a strong entropic effect arising from the energetic 
degeneracy of domain configurations and away from the 
bulk critical regime we find 



h{ 



/ 



T 

1 ln2 

2 



(9) 



which is precisely the SOS result obtained earlier by 
Duxbury and Orrick |0 and Lipowski |l^. For the di- 
agonal lattice on the other hand such entropic effects in- 
fluence the wetting but not filling phase boundary. The 
SOS limit for this lattice has not been considered before 
but is very easily constructed. Assuming the interface 
sticks to the walls (as shown) outside of the filled region 
the probability distribution for diagonal interfacial height 
follows as 



2 



P(z) = 



(10) 



where (z) ~ [h,\ — h()~^ and limr^o = 1/2 true to 
the exact Ising result. Notice that the scahng of P{z) is 
precisely the same as that predicted for the SFL regime 
of critical wetting in agreement with continuum effective 
Hamiltonian studies of filling in open wedges [Q. 

Returning to the Ising phase boundary, we note that 
the above results allow us to determine the exact R2{ip) 
even though it is not always possible to (easily) construct 
Ising lattices with arbitrary opening angle. Near the 
bulk critical temperature the Ising model recovers fluid 
isotropy and the incremental free-energy Tp may be iden- 
tified with TC0s9. Using the thermodynamic condition 
for filling (|lj) we obtain, after a little manipulation, the 
elegant expression 



i?2(V') = V2sin 



(11) 



which recovers the result for the right- angle corner 
quoted above as a special case. Note that for a corner 
with 120° opening angle we predict R2{n/3) = 1/V2 , 
which may be tested using the triangular lattice Ising 
model which naturally forms corners of this type. Also 
( |ll| ) is consistent with the general requirement for the 
open wedge {ip n/2) limit l-Rdiip) ~ (7r/2- V')^/^~"= 
which follows from (|l|) . Here denotes the critical wet- 
ting specific heat exponent which is zero for the planar 
Ising model [nSl. 
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FIG. 3. Numerical transfer matrix results for the mag- 
netisation profile in the finite-size Ising square for A*' — 24, 
T = 1.6 and different surface fields above, below and (nearly) 
at the filling phase boundary. 

To test the connection between filling and the SFL 
regime for the full Ising model we consider a finite-size 



square lattice Ising square of side length N with weak- 
ened bonds around the perimeter. The boundary spins 
are fixed to -1-1 (-1) along the bottom (top) rows and left 
(right) column respectively so that the domain wall runs 
from the upper left to lower right corner (see inset Fig. 3). 
In the limit of N tending to infinity the system decouples 
into two separate corner lattices, each of which exhibit a 
filling transition by the opposite bulk phase at the same 
temperature, Tf. For finite TV this anti-symmetric choice 
of boundary conditions frustrates the interface analogous 
to behaviour in infinite Ising strips with opposing surface 
fields d^JI^ . Far above the filling transition the interface 
stretches from the top left to bottom right corners across 
the middle of the square and has an r.m.s width of order 
^/N. In contrast, far below the transition temperature 
we anticipate pseudo phase coexistence, since the inter- 
face may be bound to either the top-right or the lower-left 
corner. The cross-over between these different types of 
interfacial behaviour occurs when Tf — T ^ ]\[-^/0o ^e- 
flecting the influence of length-scales associated with the 
filling transition. 

The conjectured connection between filling and SFL 
wetting suggests that dXT — Tf the magnetisation profile 
along the diagonal m(z) is universal and characterised 
by the same scaling function describing the profile for 
an infinite Ising strips with opposing surfaces fields at 
T — Tuj. Thus exactly a,t T = Tf and for asymptotically 
large z,N with arbitrary z/N e (0,1), we expect the 
SFL-like scaling behaviour 



m(z) — mo 1 



2z 



(12) 



where toq is the bulk spontaneous magnetisation. No- 
tice that at the filling temperature the r.m.s. width of 
the interface now scales with N compared to for 
T > Tf which reflects the additional excitations present 
near filling 0. Using numerical transfer matrix tech- 
niques 13 we have calculated the exact magnetisation 
profile in the Ising square for systems up to = 24. The 
results verify the qualitative character of the interface 
delocalisation discussed above and are in excellent quan- 
titative agreement with the scaling of the profile at Tf. 
Indeed, testing for the most linear profile turns out to be 
a highly efficient numerical method of determining the 
phase boundary up to T = 1.6 (see Fig. 2) beyond which 
the influence of the bulk critical point becomes apparent 
and larger systems sizes are required to see the scaling 
of m(z). At low temperatures the exact Ising results are 
very well described by the SOS approximation and for 
this model we have numerically computed to(z) for sys- 
tems up to A^ = 180. Again there is excellent agreement 
with the scaling prediction (11) for the profile a,t T = Tf 
and the SOS phase boundary h({T) is recovered to four 
significant figures. 

It is straightforward to derive an expression for Rd 
using mean-field (MF) theory. This serves not only as 
point of comparison with the exact Ising result but also 



3 



indicates the d dependence. We omit the details and 
only quote the final result (which is independent of d) 
obtained using Landau theory 

smij^^Ril-^) (13) 

which we anticipate to be correct for d > 4. This is 
a remarkably good approximation to the exact d = 2 
result (yielding R{tt/4) = 0.518 and R{tt/3) = 0.684) 
and strongly suggests that the dimension dependence of 
Rdi'4') is rather weak. Unfortunately no e expansion re- 
sults for the pertinent surface tensions are available to 
systematically estimate Rd in lower dimensions; though 
progress can be made using approximate non-classical 
local-functional theory which is usually a reliable treat- 
ment of bulk and surface criticality. Strictly speaking 
the theory fails for d < 3 when ip is very close to 7r/2 
since it does not account for capillary-wave like fluctua- 
tions which alter the Landau MF value of as (equal to 
zero). However, the asymptotic critical regime for wet- 
ting in three dimensions, where such effects are impor- 
tant, is extremely small and the value of i?3(V') will be 
uninfluenced by capillary- waves for all but the shallowest 
of wedges. The local functional theory we use is a sim- 
ple extension and application of the Fisk-Widom theory 
p5[ . Again we omit the details of the calculation and 
only quote the result for three dimensions JTsf which is 
conveniently written as the expansion 

, 7i?3 9Rn 

smipKi — ^(1 i + ■■) (14) 

2V6 4 160 ^ ^ ' 

valid provided i?3 is not extremely close to unity. For the 
right-angle wedge the local theory predicts i?3(7r/4) w 
0.53 which lies close to and between the d = 2 and d = 4 
(MF) results and is presumably reliable to within a few 
per cent. This may be used to estimate the location of the 
filling phase boundary for Ising model right angle wedges 
with weak surface fields which are natural generalisations 
of models used to study wetting Again, as in d = 2, 
there are two ways of constructing this wedge with a sim- 
ple cubic lattice with cross sections the same as shown 
in Fig.l (and translational invariance in the other spatial 
direction). Away from the critical region we anticipate 
that, for the analogue of the (S) lattice, entropic effects 
are important and the phase boundary is determined by 
a linear law similar to (^) . This prior knowledge of the 
phase boundary may be useful when studying the large 
scale fluctuation effects predicted for wedge filling char- 
acterised by the universal critical exponent /3o = 1/4 [|j . 

In summary we have presented exact results for cor- 
ner wetting in d — 2 and approximate results in higher 
dimensions, highlighting the universal shift of the phase 
boundary relative to wetting in the bulk critical region. 
In d = 2 numerical transfer matrix studies of the finite- 
size scaling of the magnetisation profile in Ising squares 
strongly support the conjectured relation between filling 
and the SFL regime. 
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